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Abstract 

We investigate synchronization in the presence of delay time modulation for ap- 
plication to communication. We have observed that the robust synchronization is 
established by a common delay signal and its threshold is presented using Lyapunov 
exponents analysis. The influence of the delay time modulation in chaotic oscillators 
is also discussed. 
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Synchronization phenomenon was discovered by Huygens in two coupled pen- 
dulum clocks [1]. The first description of synchronization in chaotic systems 
by Fujisaka and Yamada [2] gave new rise to scientific attention to the phe- 
nomenon and, over the past decade, synchronization in chaotic oscillators has 
been rigorously investigated not only for understanding of its role in nonlinear 
systems but for its potential applications in various fields [3,4]. Synchroniza- 
tion in chaotic oscillators is characterized by the fact that transverse variable 
X = xi— x 2 converges to zero in the limit of t — > oo, while longitudinal variable 
x = xi + X2 is still in a chaotic state. Up to now the kinds of synchronization 
[3,4,5] that have been observed include complete synchronization in identical 
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systems and phase synchronization (PS), lag synchronization, and frequency 
synchronization in slightly detuned systems, while generalized synchronization 
and generalized phase synchronization have been described in systems with 
different dynamics. 

One of the most important fields where chaos has been practically applied is 
in communication. The major concern in this field was that an encoded mes- 
sage is vulnerable to extraction by nonlinear dynamic forecasting [6], when 
it is masked by the signal from the low-dimensional chaotic system. Then it 
becomes essential to develop high-dimensional chaotic systems having a mul- 
tiple number of positive Lyapunov exponents to implement a communication 
system. In this regard, one time-delayed system received a lot of attention: 
that is, 

X = f(x(t),x(t-T )), (1) 

where r is a delay time [7,8,9,10,11]. The time-delayed system exhibits in- 
triguing characteristics with increasing tq. Despite a small number of physical 
variables, the embedding dimension and number of positive Lyapunov expo- 
nents increases as the delay time increases, and the system eventually transits 
to hyperchaos[12]. For this reason, a time-delayed system was considered to 
be the most suitable candidate for implementing communication systems. 

However, it was newly discussed the communication based on a time-delayed 
system is still vulnerable because the delay time To can be exposed by several 
measures, e.g., autocorrelation [13], filling factor[14], and one step prediction 
error, etc [14,15]. If the delay time is known, the time-delayed system collapses 
to the simple manifold in the (x,x(t),x(t — To)) space [14,15]. So the time- 
delayed system becomes quite a simple one to someone who knows the delay 
time r , and that message encoded by the chaotic signal can be extracted by 
the common attack methods [6]. 

Accordingly, the most crucial questions to address regarding to the application 
of time-delayed system to communications come to be "Is there any scheme 
which hides the delay time from an eavesdropper", and "Is it possible to find 
out a robust synchronization regime under the scheme". In this paper, we con- 
sider a generalization of time-delayed systems in which the delay time is not 
constant but modulated in time, which was introduced to study the general 
property of the delay differential system [16]. Applying the generalization to 
two coupled chaotic systems, we focus our discussion on answering the two 
questions above. Thus, we would show that there exists a robust synchroniza- 
tion regime even if the delay time is modulated and that the imprints of delay 
time in time series are completely wiped out by the modulation of delay time. 

First, to explain the delay time modulation (DTM), we consider a single lo- 
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Fig. 1. Autocorrelation functions and filling factors at A = 4.0, (3 = 0.7: (a) 
and (c) constant delay time tq = 30; (b) and (d) delay time modulation with 
g(t, To) = (0.95to/2) sin(i) + r /2 and tq = 60. The primary peaks in (a) and (c) 
show the constant delay time tq. 



gistic map with time-delayed feedback, and observe what an effect the mod- 
ulated delay time r has on the characteristics of the time series: x(t + 1) = 
\x(t)(l-x(t)), where x(t) = (1 - P)x(t) + (3x{t - r) with r(t) = g(t,r ). The 
modulation function of the delay time is given by: 



g(t,T ). 



(2) 



If we take g(t, To) as a constant, the system returns to the usual time-delayed 
system. As we shall see, the nontrivial choice of the modulation function not 
only has a crucial effect on the characteristics of synchronization, but also 
remarkably transforms the characteristics of the time series. Figure 1 shows 
autocorrelations and filling factors [17] for two different modulations. In the 
case of the constant delay time (Fig. 1 (a) and (c)), one can easily notice that 
the peaks indicated by arrows exactly correspond to the multiples of the delay 
time To- On the contrary, when the modulation is sinusoidal (Fig. 1(b) and 
(d)) the imprints of delay time are completely wiped out, and the detection of 
delay time seems to be impossible. We emphasize that the disappearance of 
the imprints of delay time is caused by the fact that the delay time is changing 
in time. Accordingly, Figure 1 clearly shows that the results do not depend on 
the method of the detection. 



To analyze synchronous behaviors, we consider a system of two coupled Lorenz 
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Fig. 2. Projective trajectories of Eq. (3) when j3 = 0.93 and to = 0.005. (a) in (x, y) 
space (b) in (y, z) space. 

oscillators with common time-delayed feedback: 



xi = a{yi - xi), 

yi = -XiZi + rxi - y u 

zi = xiUi ~ bz u 

r = g(t,T ), 

%2 = cr(y 2 -x 2 ), 

y 2 = -x 2 z 2 + rx 2 - y 2 , 

z 2 = x 2 y 2 - bz 2 , 



> Transmitter 



(3) 



> Receiver 



(4) 



where a = 10.0, r = 28.0, and 6 = 8/3 and x ± = (1 - P)x 1 (t) + Py^t - r) 
and x 2 = (1 — f3)x 2 {t) + (3yi(t — r). Here we take the modulation function 
as g(t, t ) = °' 9 | r ° sin(cjt) + y and r = 60. We note that the systems are 
coupled by common signal (3yi(t — r). Figure 2 shows the trajectories in the 
phase space with time-delayed feedback. One can see that the shape of the 
attractor is so seriously deformed that it is impossible to identify the original 
shape of the Lorenz attractor. This deformation of the attractor is the typical 
characteristic of the time-delayed system which stems from the transition to 
hyperchaos by time-delayed feedbacks [7,14]. 

The synchronization transition can be easily detected by observing the relative 
motion of the two oscillators. The equation of motion describing the relative 
motion is obtained from Eqs. (3)- (4): 



X = a(Y-(l-P)X), 



Y = —(1 — P)(xiZ + Xz\ — XZ) 
-P yi (t-r)Z-r(l-P)X-Y : 



(5) 
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Fig. 3. Two largest TLEs as a function of the coupling (3 and the modulation fre- 
quency uj. We have added the uniform random noise of amplitude |£| ~ 10 -8 only 
to x variable of Eq. (3) in order to remove abrupt synchronization regime: (a) the 
largest exponent, where the border line between gray and dark gray region on ((3, to) 
plane corresponds to the transition threshold to CS; (b) the second largest exponent, 
where the border line corresponds to the transition threshold to PS. 




Fig. 4. Contour plot of synchronization regimes in (u,(3) space. A and B are the 
reference points to show the time series. CS: Complete Synchronization, PS: Phase 
Synchronization, AD: Amplitude Death. 



Z = (l-/3)(x 1 Y + Xy 1 -XY) 
+P yi (t-r)Z-bZ, 

where X = x\ — x 2 , Y = y\ — y 2 , and Z = z\ — z 2 . The above equations are 
nonautonomous in themselves and, so, by iterating the above equation with 
Eqs. (3)-(4), we obtain the transverse Lyapunov exponents (TLEs) [18] which 
enable us to determine the synchronization threshold. 
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Fig. 5. Time series of x\ — X2 and the delay time r as a function of time: (a) point 
A in Fig. 4 ((3 = 0.87, u = 0.005). (c) point B in Fig. 4 ((3 = 0.93, u = 0.005). (c) 
modulated delay time when to = 0.005. 



Figure 3 shows two largest TLEs as a function of (3 and u. One can clearly 
see the threshold to CS on ((3, uj) plane of Fig. 3 (a) where the largest TLE 
becomes negative. Also the threshold to PS, where the vanishing TLE be- 
comes negative, is plotted on ((3,ui) plane of Fig. 3 (b). The synchronization 
regimes are presented in Fig. 4 as a contour plot. Contrary to a common intu- 
ition, the CS regime becomes wider as the modulation frequency u increases, 
even though the regime is smeared with an amplitude death (AD) [19,20]. 
The amplitude death phenomenon due to DTM in our study is similar to the 
transition to periodic orbit in Ref. [16] in the sence that in both cases the 
systems lose their exponential instability by the interaction. However differ- 
ent apect is that our system converges to the fixed point without oscillation 
instead of the periodic orbit. These effects show that DTM can play a role 
of stabilizing scheme such as chaos control[10]. The appearance of the wider 
synchronization regime in Fig. 4 implies that the scheme will be effective to im- 
plement real communication systems, since the higher modulation frequency 
widens the synchronization regime and it would enable us to construct more 
of communication channels. 



Figure 5 shows the temporal behaviors of the difference between x\ and x 2 and 
the modulated delay time r. The temporal behaviors below (point A in Fig. 4) 
and above (point B in Fig. 4) the threshold are presented in Fig. 5 (a) and Fig. 
5 (b), respectively. They clearly show the threshold behaviors of CS at two 
points. The delay time r is modulated as shown in Fig. 5 (c). It is interesting 
to see the autocorrelations and filling factor for different modulations. For the 
constant delay time, the delay time is exactly indicated by the position of the 
dips pointed by arrows in Fig. 6 (a) and (c). In the case of DTM (Fig. 6 (b) 
and (d)), the dips disappear. 

In general introducing a d- dimensional chaotic oscillator, we can consider a 
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Fig. 6. Disappearance of the dips in autocorrelations and filling factors at (5 = 0.93: 
(a) and (c) constant delay time, tq = 30; (b) and (d) modulated delay time 
(oj = 0.005). 

chaotic modulation of delay time: r = g(x, To) where x = -F(x) and e is 
a time scaling parameter to control the rate of modulation like u in sinu- 
soidal modulation. We can expect that the chaotic modulation will not work 
a fundamental change in contrast with the sinusoidal modulation, because 
the delayed feedback x(t — r) hardly correlated with x(t) already even in the 
sinusoidal modulation as shown in Fig. 6 (b). However, in the application 
to communication, the chaotic modulation would enhance the security of the 
system because it may increase the complexity of the attractor. It is worth 
discussing the physical implementation of the systems described by Eq. (3)- 
(4). We think that such a situation can be observed in optical systems. For 
example, one can consider two laser systems coupled by reflected laser beam 
from a vibrating mirror. In that case, Lorenz equation describes a laser system 
and the vibrating mirror plays a role of modulation function. 

In conclusion, we have studied the effects of DTM in chaotic oscillators. We 
have observed that the imprints due to the constant time delay are completely 
wiped out by the modulation in the chaotic map and flow. We have also demon- 
strated that the robust synchronization can be established in the presence of 
DTM. We expect this effects to be useful in constructing real communication 
systems. 

We thank M. S. Kurdoglyan for helpful discussions. This work is supported 
by Creative Research Initiatives of the Korean Ministry of Science and Tech- 
nology. 
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